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Abstract. We review a few rigorous and partly unpublished results on the regularisation 
of the stress-energy in quantum field theory on curved spacetimes: 1) the symmetry of the 
Hadamard/Seeley-DeWitt coefficients in smooth Riemannian and Lorentzian spacetimes 2) the 
equivalence of the local (^-function and the Hadamard-point-splitting procedure in smooth static 
J> I spacetimes 3) the equivalence of the DeWitt-Schwinger- and the Hadamard-point-splitting pro- 

■ cedure in smooth Riemannian and Lorentzian spacetimes. 

■ 1 Introduction 

(N 

^ I The "semiclassical" Einstein equation 



G^, = 87rG{:T^,:)n, (1) 

is obtained from the "classical" Einstein equation by replacing the classical stress-energy 
tensor T^^, with the expectation value of the quantum stress energy tensor : T^j, : in a 
quantum state Q. There are several angles from which one can approach the question 
whether this equation is sensible and well-defined. On the one hand, one can argue that 
it can be obtained from a theory of (perturbative) quantum gravity in a suitable limit 
where gravity becomes classical, but matter is still quantized, hence the name; see e.g. 
|FlWa96j for a review of this issue. On the other hand, if we combine our understanding 
of e.g. radiation being described by quantum field theory on the fundamental level and 
the fact that radiation does have an influence on spacetime curvature for all we know 
from cosmological observations, we can safely claim that the latter observations indicate 
that (something like) the semiclassical Einstein equation must hold at least in a suitable 
regime, although it is certainly not a fundamental equation, as it equates a deterministic 
classical quantity with a probabilistic quatum expression. 

Even if one accepts one or both of the aforementioned physical motivations to consider 
the semiclassical Einstein equation, questions about the well-posedness of this equation 
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remain. We usually assume that spacetime is regular outside of black hole or big bang sin- 
gularities, so (: T^u ■)n should better be regular as well. However, as the classical expression 
for T^j, is (at least) quadratic in the field, one can expect that the naive quantum stress 
energy tensor is divergent and has to be regularised in some way. Indeed, various schemes 
to regularise the quantum stress-energy tensor have been proposed, including some more 
familiar from flat spacetime like dimensional regularisation |Bro76] . (^-function regularisa- 
tion |DoCr76[ IHaw77] or Paul-Villars regularisation |BiDa82] and also ones which are less 
familiar in that context like point-splitting regularisation |Chr76t [Chr78t rWal77t IWal78a] . 
Naturally, the question arose which regularisation scheme was the "correct" one, and a 
satisfactory answer was first given in |Wal77] (and later generalised in [HoWa05] ) . In 
|Wal77] . Wald proposed a minimal set of well- motivated conditions which any sensible 
regularisation scheme should satisfy. We state these conditions, the "Wald axioms" , here 
in their modern form |Wal95l IBFV03l[HoWa05] . 

1. The commutator of the regularised stress-energy tensor : T^^{x) : with any product of 
fields (j){xi) ■ ■ ■ 4){xn) equals their commutator with the non- regularised stress-energy 
tensor T^i,{x). 

2. ■.T^y{x) : transforms covariantly under diffeomorphisms and does not depend on the 
metric and its derivatives ai y ^ x. 

3. (The expectation value of) -.T^^^x): has vanishing covariant divergence. 

4. In Minkowski spacetime and in the Minkowski vacuum state the expectation value 
of -.T^^i^x): vanishes. 

5. (The expectation value of) : T^^{x) : does not contain derivatives of the metric of 
order higher than two. 

While the first condition makes sure that the regularisation scheme is "minimal" and 
"state-independent" , as the divergence in the non-regularised stress-energy tensor is pro- 
portional to the identity, the second assures that is is local and purely geometric, and 
the third follows from the fact that G^y has vanishing covariant divergence as well. The 
fourth condition, introduced to guarantee that the regularisation prescription coincides 
with "normal ordering" in flat spacetime, has to be omitted if one wants to have a non- 
vanishing cosmological constant (on the right hand side of the Einstein equations) |Ful89] . 
Finally, the last condition was meant to avoid the appearance of unstable, "runaway" so- 
lutions to the semiclassical Einstein equation. Although it turned out that this condition 
can not be fulfilled and that such unstable solutions can not be avoided, one can still 
analyse the subclass of stable solutions in a controlled manner in special cases, see e.g. 
[Statol IVil85l IDFP081 iKokOOl [PmTOl [PHMPlOj . 

Wald has conjectured already in [Wa79j that C-function regularisation, dimensional 
regularisation and the various point-splitting regularisation schemes satisfy the essential 
axioms 1., 2. and 3. and thus all have a raison d'etre although they give slightly different 
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results. In fact, one can infer directly from the axioms 1., 2. and 3. that their solution 
can not be unique, but is subject to a finite regularisation freedom parametrised by 4 
divergence-free tensors constructed out of the metric and the Riemann tensor |Wal77t 
IHoWaOS] : in the absence of a fundamental theory of quantum gravity and matter, this 
freedom can only be fixed by comparison with experiments [DFPOSl IDHMPIO] . So Wald 
could not only motivate why most regularisation prescriptions on the market should be 
legitimate but also explain why they naturally give different results. Nevertheless, the 
understanding of the renormalisation ambiguity still seems to be a topic of active research, 
see e.g. |ABPS06j . 

In this work we shall be concerned with proving parts of Wald's conjecture, namely, 
the rigorous relation between three different regularisation schemes, obtained in partially 
already published |Mor98t IMor99t IMorOO] and partially yet unpublished |HaclO] previous 
works of the authors. On the one hand, we shall review the equivalence between the local 
^-function regularisation and the Hadamard point-splitting regularisation in Riemannian 
( "Euclidean" ) manifolds and the compatibility of the local ^-function scheme with Wick- 
rotation whenever applicable. This scheme is only meaningful in Riemannian manifolds 
because it is based on spectral properties of the field equations, which strongly depend on 
the metric signature. Hence, a rigorous comparison with the Hadamard scheme is only 
possible in the Riemannian context a priori and can be taken over to the Lorentzian setting 
by means of a Wick-rotation only in static Lorentzian spacetimes. The situation is com- 
parable in the second case we shall discuss, the well-posedness of the DeWitt-Schwinger 
prescription. To wit, like the local ^-function scheme and for the same reasons, the 
DeWitt-Schwinger prescription is only well-defined in Riemannian manifolds. Neverthe- 
less, the computations of the regularised stress-tensor for free quantum fields of arbitrary 
spin in Lorentzian spacetimes in |Chr78t IChDu79] . obtained with the DeWitt-Schwinger 
prescription, are now part of the standard literature |BiDa82j . although it was clear that 
their derivation had not been rigorous in Lorentzian spacetimes and that their legitimation 
via a Wick-rotation was not possible in general Lorentzian spacetimes (see for instance 
the comments at the beginning of section 6.6 in |BiDa82] ). Moreover, as we shall explain 
in the main body of this work, the DeWitt-Schwinger prescription to regularise the ex- 
pectation value of the stress-energy tensor fails to take into account the state-dependence 
of this object and thus only gives the "geometric piece" of (: T^,^ this, however, is 
sufficient to determine the so-called conformal anomaly |Duf77| IWal77] . Hence, although 
the results of |Chr78t IChDu79] on the conformal anomaly have been confirmed with rigor- 
ous calculations for Klein-Gordon |Wal77t IMor03] and Dirac |DHP09] fields by now, it is 
still important to know whether the DeWitt-Schwinger scheme is generally legitimate on 
arbitrary Lorentzian spacetimes and if one can formulate it in a way which does take into 
account the full state dependence of (:T^^:)a; we shall prove both of these aspects in what 
follows. Although our proof is based on the direct relation between the so-called Seeley- 
DeWitt coefficients and the Hadamard coefficients, the DeWitt-Schwinger prescription 
still differs from the Hadamard scheme in that it does not consistent of just applying a 
suitable differential operater to the regularised two-point function. Our results imply that 
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the local ^-function and the DeWitt-Schwinger prescription both satisfy Wald's axioms 
just like the Hadamard regularisation scheme, though, in the (^-function case, only on 
Lorentzian manifolds where a Wick-rotation is possible. 

As the Seeley-DeWitt/Hadamard coefficients are of such an utmost importance in the 
point-splitting regularisation prescriptions (and in non-commutative geometry |CoCha97] ). 
we complement our analysis of the different regularisation schemes by reviewing a proof 
of their symmetry. Although this symmetry had been heavily used and expected to hold 
because these coefficients also appear in the symmetric Green's functions of the field 
equations |Fri75t IGar64] . a rigorous proof was missing since the coefficients appear only 
in infinite series whose convergence is unclear except in analytic spacetimes. The proof 
of |Mor99t IMorOOj proceeds in three steps. First the symmetry is proven in analytic Rie- 
mannian manifolds by employing the symmetry of the field equations. Then, the proof 
is extended to analytic Lorentzian spacetimes by means of a local Wick rotation, i.e. one 
"rotates" the metric rather than the coordinates. Finally, the symmetry is established in 
general non-analytic Lorentzian spacetimes by locally approximating in a controlled way 
smooth metrics with analytic ones. 

Our paper is organised as follows. After introducing a few preliminary notions in 
section 2, we continue in section 3 by discussing the expansion of the heat kernel in terms 
of the Seeley-DeWitt coefficients and the expansion of the two-point function of quantum 
field in terms of the Hadamard coefficients, followed by a proof of their symmetry in 
section 4. In sections 5 and 6 we review the proof of the equivalence of the local (^-function 
and the Hadamard regularisation scheme in Riemannian manifolds and static Lorentzian 
manifolds and the well-posedness of the DeWitt-Schwinger regularisation prescription in 
Lorentzian spacetimes, respectively. The paper ends with our conclusions in section 7. 

For simplicity and to increase the accessibility of this work, we consider only the real 
Klein-Gordon field in four dimensions. However, some of the results we review have been 
proven in arbitrary dimensions |Mor98t IMor99t IMorOOj , and we have no reason to doubt 
that they can be extended to fields of higher spin as well. 

In the following a manifold (including bundles) is supposed to be smooth (i.e. C°°), 
unless specified otherwise, Hausdorff and second countable (to assure paracompactness). 
Vector fields, tensor fields, metrics and sections are supposed to be smooth as well unless 
specified otherwise. We choose the signature (— , -|-, -|-, -|-) for Lorentzian metrics. 

A spacetime is a time oriented four- dimensional manifold equipped with a smooth 
Lorentzian metric. In the text, Riemannian and Euclidean are synonyms. 

2 Gaussian states, static spacetimes, Wick rotation 

Throughout we suppose to consider globally hyperbolic spacetimes (M, g) only and we 
moreover stick to the case of a real scalar field (p propagating in {M,g), satisfying: 

P0 = , P := -V^V^ + V = -D + V , (2) 
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the covariant derivative V being always associated with the metric g and V being a 
smooth scalar function of the form 

V{x) := ^R + m^ + V'{x) . (3) 

Above ^ is a real constant, R is the scalar curvature and m > the mass of the particles 
associated to the field. The domain of P is the space of real-valued C°° functions with 
compactly supported Cauchy data. A quasifree state - also known as Gaussian state - Q 
is unambiguously defined by choosing a Wightman two-point function {(j){x)(f){y))n. The 
GNS theorem or the standard Wick expansion allows to construct a corresponding Fock 
space realization of the theory. 

A globally hyperbolic spacetime (M, g) is said to be static if it admits a time-like 
Killing vector field X := dt normal to a smooth spacelike Cauchy surface S. Consequently, 
there are (local) coordinate frames {x^ , , x"^ , x^) = (t,x) where goi = {i = 1,2,3) and 
dtg^iu = and x are local coordinates on E. When discussing the static case, we also 
assume that V satisfies dtV = so that the space of solutions of ([2]) is invariant under 
t-displacements. 

In a globally hyperbolic static spacetime it is possible to construct a t-invariant Gaus- 
sian state fl by performing the Wick rotation, i.e. obtaining the Euclidean formulation 
of the same QFT. This means that one can pass from the Lorentzian manifold (M, g) 
to an associated Riemannian manifold (ME,g^^^) by the analytic continuation t ir 
where r G M. In this way, dt gives rise to another Killing vector dr in Me- Me can 
alternatively be defined by assuming that the orbits of the Euclidean time r are closed 
with period (3, including the case r G M as /3 = -|-oo. Within this approach |FuRu87] . the 
two-point function of Q is completely determined by a proper Green's function (in the 
spectral theory sense) S{t — t',x,x') of a corresponding self-adjoint extension Pe of the 
operator 

P'^ := _^iE)^iE)a ^ _ _^{E) ^ y^^^ . C^^Me) ^ L\Me, dg,E)X) . (4) 

S{t — t', X, x') is periodic with period (3 G (0, +00) U {+cxd} in the t — r' entry and it is 
called the Schwinger function of fi. S is the integral kernel of P^^ when Me is compact 
and P^ > |Wa79j and, in that case, Pe = P'e, that is Pe is essentially self-adjoint. In 
general a t-invariant state Q is determined by the choice of a self-adjoint extension of the 
Wick rotated Klein-Gordon operator P'^. When < f3 < +00, the value T = 1/(3 has to 
be interpreted as the temperature of a state fl because the Wightman two-point function 
{(j){x)(f){y))fi satisfies the KMS condition at the inverse temperature (3. 

3 Hadamard / Heat-kernel Seeley-DeWitt coefRcients 

In order to introduce the Hadamard and Seeley-DeWitt coefficients we need a few defini- 
tions and results from bitensor calculus. 

If VM and WM are vector bundles over a manifold M with typical fibers constituted 
by the vector spaces V and W respectively, then we denote by VM Kl WM the exterior 
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tensor product of VM and WM. VM Kl WM is defined as the vector bundle over M x M 
with typical fibre V ® W. The more familiar notion of the tensor product bundle VM ® 
WM is obtained by considering the pull-back bundle of VM Kl WM with respect to the 
map M 9 X (-> (x, x) G M^. Typical exterior product bundles are for instance the tangent 
bundles of Cartesian products of M, e.g. T*MMT*M = T*M'^. A section of VMMWM 
is called a bitensor. We introduce the Synge bracket notation for the coinciding point 
limits of a bitensor. Namely, let i? be a smooth section of VM Kl WM. We define 

[B{x,y)] := \imB{x,y) . 

With this definition, [B{x,y)] is a section of VM ® WM. In the following, we shall 
denote by unprimed indices tensorial quantities at x, while primed indices denote tensorial 
quantities at y. Moreover, covariant derivates shall be denoted by the usual abbreviated 
notation 

B;^ := . 

As an example, let us state the well-known Synge rule, proved for instance in |Chr76] . 

Lemma 1 Let B be an arbitrary smooth bitensor. Its covariant derivatives at x and y 
are related by Synge's rule. Namely, 

[B,^,] = [B],,-[B,,]. 

Particularly, let VM be a vector bundle, let fa be a local frame ofVM defined on O <Z M 
and let x, y E O. If B is symmetric, i.e. the coefficients Bab'{x,y) of 

B{x,y):=B-'\x,y) fa{x)®fy{y) 

fulfil 

B^''\x,y) = B''\y,x), 

then 

[B.,,,] = [B,,] = \[B],,. 
Other relevant properties of [■] can be found in e.g. |Ful89] . 

3. 1 Hadamard coefficients, Hadamard states and point-splitting regularisation 

An important tool of the Hadamard point-splitting renormalisation prescription in curved 
spacetime are the so-called Hadamard coefficients. They appear in the explicit expression 
of the Hadamard parametrix. That notion, in turn, has its physical reason in the definition 
of a relevant class of quantum states called (locally) Hadamard states |KaWa91] IRadOGaj 
IRad96bt IWal95] . To discuss them, we consider a Gaussian state Q for a real scalar field 
on a globally hyperbolic spacetime (M, g) as before. 
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Definition 2 We say that {(f){x)(f){y))n is of local Hadamard form, andVt is locally 
Hadamard, if: 

•= li"^ ^ y) + w{x, y)) , 
where x, y are taken in any geodesically convex neighborhood. Above: 

a,{x, y) ■= (t(x, y) + 2ie{T{x) - T{y)) + (5) 

where T is any local time coordinate increasing toward the future, and (j{x,y) denotes 
one half of the signed squared geodesic distance of x and y (well-defined in a geodesically 
convex neighborhood) and X > a renormalization length scale. Finally w is a smooth 
biscalar determined by the state up to re-definition of X. 

The bidistribution hi, is the Hadamard parametrix and the coefficients (bi-scalars) 
Vn are called Hadamard coefficients. 

The convergence of the series in the expression of is just asymptotic. However a proper 
point-wise convergence can be obtained by introducing a set of suitable cutoff functions 
(see |Fri75t IHoWaOl] for details). The parametrix /i^ solves the Klein-Gordon equation 
in each argument separately up to smooth terms. From that one obtains the following 
recursive differential equations (actually they do not depend on the signature of the met- 
ric) already studied by Riesz in his seminal paper |Rie53] (see also |Gar64] ) . determining 
u = u{x, y) and all the Hadamard coefficients f„ = Un(a^, y) in a fixed geodesically convex 
neighborhood when assuming [u] = 1. 

2m.^(t.^ + (D^a - 4)m = , (6) 
-P^u + 2vo-f,cr.^ + {D^cr - 2)vo = , (7) 
-P,Vn + 2{n + l)vn+i-f,(T.^+{{n + l)n,a + 2n{n + l)){vn+i) = Vn > (8) 

In particular, it turns out that u coincides with the square root of the so called VanVleck- 
Morette determinant |Ful89l IMor99] . As a consequence, the following relations can be 
proven (see e.g. [Mor03j ) : 

Lemma 3 The following identities hold for the Hadamard parametrix hi;{x, y) 

[P^h,] = [Pyh,] = -6[v,] , [iP.h)J = [iPyh,)y] = -A[v,].^ , 

[{P.K)-y] = mA)J = -2[v^],,. 

These identities have a great deal of effect concerning the point- splitting renormalisation 
of the stress energy tensor evaluated on a Hadamard state Q. It is essentially defined by 
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subtracting the universal Hadamard singularity from the two-point function of il, before 
computing the relevant derivatives. 

{■Tf,u - (x))^ = [D(^.^)^^{x,y) {{(f){x)^{y))^ - hs{x,y))] (9) 

Di_i.„{x, y) is the following second order partial differential operatoi0 (cf. |Mor03j eq. (10), 
and |HaclOj . |ElGoll] where some minor misprints have been corrected), 

D^,{x,y) := Dl^J'{x,y) + ^g^,P, 
D^P{x, y) := (1 - 20^?r' V^V.. - 2eV^V. + ^G,, 

+ g,, |2en. + {2^ - <'V^Vp, - ^m^l 

Here covariant derivatives with primed indices indicate covariant derivatives w.r.t. y, g^' 
denotes the parallel transport of vectors along the unique geodesic connecting x and y, 
the metric g^i, and the Einstein tensor G^^ are considered to be evaluated at x and we 
have assumed V = for simplicity. 

The form of the "canonical" piece follows from the definition of the classical 
stress-energy tensor, while the last term ^g/^uPx has been introduced in |Mor03] and it 
gives no contribution classically, just in view of the very Klein-Gordon equation satisfied 
by the fields. However, in the quantum realm, its presence has important consequences 
due to Lemma |3] because the Hadamard parametrix satisfies the Klein-Gordon equation 
only up to smooth terms and thus the above definition without this additional term would 
not yield a conserved stress-tensor expectation value. As a matter of fact the following 
theorem is valid ( |Mord3] Theorem 2.1). 

Theorem 4 The stress-energy tensor expectation value defined by ^ fulfils the Wald 
axioms 1., 2., 3. for every value of X > and, with a suitable choice for the scale A > 0, 
it satisfies 4- For all A it displays the trace anomaly when the conformal coupling ^ = 1/6 
and m = are chosen: 

0^" i-T ■ ix)) = 

To conclude we remark an important theoretical result: if a globally hyperbolic space- 
time is static and Q is defined as in sec. [2] through the Wick rotation procedure, fl turns 
out to be Hadamard |Wal95| IMor03] for < f3 < +00. In particular the Minkowski 



vacuum, that can be constructed in that way, is Hadamard. 



3.2 Seeley-DeWitt coefficients 

We now introduce another class of biscalars called the Seeley-De Witt coefficients. They 
appear if on tries to regularise both the Euclidean and Lorentzian one-loop effective action 

^Often, the symmetrised version of is considered. However, since {4>{x)4>{y)) — hg{x, y) is already 
symmetric, it is not necessary to symmetrise the differential operator. 
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by means of the asymptotic expansion of the heat kernel. Actually the notion of heat kernel 
is properly defined only in the Euclidean case, though the terminology has been adopted 
for the Lorentzian case as well. For the moment we focus on the mathematical issues only, 
and just from a very formal point of view, while we leave aside the physical motivations 
leading to the Euclidean approach (especially the definition of the considered quantum 
state) as well as precise mathematical definitions that will be discussed in a subsequent 
section. We just say that a formal motivation for passing from the Lorentzian to Euclidean 
side and vice versa relies upon the natural extension of the Wick rotation procedure 
discussed in sec. [2] when dealing with static spacetimes. On a Euclidean manifold {M,g), 
the Euclidean effective action S^a is defined via an Euclidean path integral as 



A being the Laplace-Beltrami operator associated with g and V" : M — M being some 
given smooth function such that P > 0, is the Euclidean Klein-Gordon operator defined 
on (M, (?). [del)] is supposed to denote some (formal) measure on the space of field config- 
urations. The argument of the exponential in the integrand is the classical (Euclidean) 
Klein-Gordon action taking boundary terms into account. On Lorentzian manifolds, for- 
mally speaking, one would put an imaginary unit in front of both actions in the above 
formula. In that case, however, the rigorous heat-kernel theory cannot be applied since it 
strongly relies upon the fact that A is (strongly) elliptic and several regularity results (in- 
cluding convergences of expansions) generally cease to hold true. See |Wa79t [Ful89t IMor99j 
for comments on these issues. 

The relevance of the effective action in our case is the fact that one can formally define 
the expectation value of the regularised stress-energy tensor in analogy to the classical 
case as 



The natural question which arises is where the state (fi) dependence of the right-hand 
side is hidden. Indeed, one can interpret the situation in such a way that the measure 
is state-dependent |Mor03] . We shall consider these issues in forthcoming sections. 
To compute (: T^y in this picture, one therefore needs to compute the mentioned 
functional derivative with respect to the metric. As we are going to see, already Ses is 
affected by divergences and thus a regularization process is necessary. If the regularisation 
is such that the effective action is diffeomorphism-invariant, than the resulting functional 
derivative will automatically lead to a conserved {:T^u:)n. 

Let us proceed to understand why the effective action is divergent and how to regularise 
it. Evaluating the above integral as if we were in a finite dimensional vector space, thus 





(10) 




(11) 
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computing a Gaussian integral, one finds (up to a constant) 



e-^^ff = (detP)-^ =e-^^'~'°sP^ (12) 

and, hence, 

5'eff = ^ / [^ogP]{x,x) := j dgX Lcff(x) , 

M M 

where the indentity (vahd in finite dimensional cases) log det P = Tr log P has been 
used, the trace is evaluated by integrating the kernel of the operator log P at coinciding 
points, and L^s is interpreted as the effective Lagrangean. The above integral is certainly 
diverging, as one can already infer from the identity (that holds true for M compact 
when replacing P in ( ITOj) for its unique self-adjoint extension P and referring to the weak 
operator topology) 

p 

log P = lim ^ — [ ds 1- (7 — log e)I 

where 7 denotes the Euler-Mascheroni constant and I is the identity operator. We discard 
the divergent term proportional to I, that is trace class only in finite dimensional spaces, 
and define (where again, to be rigorous, one should replace P by P) 

Leff(a;) := - / ds ^ . (13) 

The integral kernel [e~'^^]{s,x,y) of the appearing exponential of P is called the heat 
kernel because it satisfies the heat equation 

(ds + P.) [e-^T {s,x,y) = 0. 

It is possible to prove that, at least when M is compact and P > or, more weakly, it 
is bounded below (but there are several other generalizations, see |Wa79t IMor99] also for 
references) the heath kernel exists, is an integrable function and, in geodesically convex 
neighborhoods, the following expansion is valid: 

r PI -1- ^{^jV) V — ^ 

[e~* \{s,x,y) = j^ — e ~ 2_^(^n{x,y)s^ + 'smooth biscalar' (14) 



n=0 

00 



1 a(x.y) 2 V — ^ 

e 2^ ^ y^an{x,y)s"' + 'smooth biscalar' 



^ ' n=0 

The convergence of the series is "asymptotic" in a certain sense discussed in sec. 1.2 
of [Mor99] (and it is not asymptotic in the Lorentzian case). The functions a„(x, y) are 
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smooth biscalars completely determined by the recursive differential equations |Ful89 
IMor99j 



ao = u, 



P^ttn + cr^a„+i;^ + -Aa + n - 1 a„+i 







(15) 



with u, again, denoting the square root of the VanVleck-Morette determinant and a 
being related to a„ via 



In the Lorentzian framework, formally speaking, one gets the same expansion and the same 
coefficients, though defined with the Lorentzian metric and, in particular A is replaced 
by □ in the recursive equations above. 

The reason why we have displayed the expansion (HM in two versions is the following: 
the version in terms of a„ is more convenient in the regularisation procedure, because 
the appearing avoids potential infrared singularities in the integral with respect to s 
present in L^s- In contrast, the version given in terms of a„ is important to show the 
relation between the Hadamard coefficients Vn and the a„. Namely, a short computation 
and comparison with the scalar Hadamard recursion relations discussed in subsection 13.11 
reveals the well-known identity 



4 The symmetry of the Seeley-DeWitt/Hadamard coefRcients 

An intriguing issue concerns the symmetry of the above introduced functions a„(x, y) (or, 
equivalently, f„(x, y)) under interchange of their arguments x and y when (M, g) is smooth 
it and either Riemannian or Lorentzian. The symmetry is by no means evident in view 
of the fact that their definition, relying on (|8]), ([7]) and ([6]) is highly non-symmetric when 
swapping x and y because derivatives act on x only. In |Wal78a] . such a symmetry was 
(indirectly) argued to hold for the analytic case. In |FSW78] . the symmetry was argued to 
hold for the C°° case case and used to simplify some technical computations. Nevertheless, 
these papers did not report the corresponding proof. In the literature concerning the 
point-splitting procedure successive to |Wal78a] . e.g. |B086] . this symmetry is assumed 
to hold implicitly. The issue was analysed in |Mor99l IMorOO] and a proof was presented 
both for the Euclidean and the Lorentzian case. In the following we review the main ideas 
of the proof, and refer to |Mor99l IMorOO] for technical details. 

Theorem 5 The Seeley-DeWitt/Hadamard coefficients a„(a;, y) defined in a geodesically 
convex neighbrhood of a manifold {M, g) , either Riemannian or Lorentzian, associated 
with the Klein-Gordon (resp. Riemannian or Lorentzian) operator P (with V smooth but 
arbitrary) are symmetric if swapping the arguments x and y. 




j=0 



(16) 




(17) 
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Sketch of proof. First of all, given a Riemannian manifold M and a geodesically convex 
neighborhood O C M, where the are defined, one modifies M and V outside O, in 
order to produce a compact manifold M' D O quipped with an essentially self-adjoint 
operator P' bounded below and coinciding with P in O. This does not affect the defini- 
tion of the ttn since they depend on the local geometry in O only. As e~'^^' is self-adjoint, 
it follows that [e~'^^']{s,x,y) = [e~'^^']{s,y,x) for x,y E O. The convergence properties 
of the expansion f|T^ |Mor99j assure in this case that {an{x,y) — an{y,x)) \x=y = 
for any multiindices referring to any fixed coordinate frame on O. If is a real 

analytic function of the said coordinates, the coefficients an{x,y) turn out to be analytic 
functions of the coordinates on O x O and, just due to the Taylor expansion, one has 
an{x,y) = an{y,x) in a neighborhood of the diagonal of the open set C x O and thus 
everywhere thereon since it is connected. Whenever the smooth metric g is not real ana- 
lytic, it is however possible to construct a class of analytic metrics gt, depending on the 
parameter t G (0, 1), uniformly converging to g as t ^ 0. In view of the structure of the 
equations (|H]), dZ]) and (jS]), exploiting a suitable version of the theorem on the analytic de- 
pendence of parameters for the solutions of a class of differential equations, one can show 
|Mor99] that the associated coefficients an\x,y) satisfy as expected an\x,y) — )■ an{x,y) 
as t — >■ 0, assuring the symmetry of an{x,y) since an\x,y) = an\y,x) if t > 0. 

Let us pass to the case where (M, g) is Lorentzian. It is now enough to establish 
the symmetry property for the analytic case because the extension to the smooth case 
is identical to that in the Euclidean case. So we assume that in a geodesically convex 
neighborhood O G M the metric g and thus, the coefficients an{x,y), are real analytic 
functions in a suitable coordinate frame thereon. As proved in |MorOO] . it is possible to 
think of (9 as a real section of a complex manifold U equipped with a complex analytic 
metric h whose restriction to O is nothing but g. The pivotal fact is that this complex 
manifold admits another real section Oe whose metric gE obtained by restriction of h is 
Euclidean and Oe is geodesically convex as well. In this sense the real analytic Lorentzian 
metric g is the complex analytic continuation of the real analytic Euclidean metric gE- 
This procedure realizes a local Wick rotation, depending, however, on the chosen coor- 
dinate frame. Once again, exploiting a suitable version of the theorem on the analytic 
dependence of parameters for the solutions of a class of differential equations (on complex 
manifolds now), one verifies that the function an{x, y) — an{y, x) in O is the analytic con- 
tinuation of a corresponding function alf\xE,yE) — cih\yE,XE) in O. Since the latter 
is the zero function, as previously established, we have an{x,y) = an{y,x) in O for the 
Lorentzian (analytic) case, too. 

5 The equivalence of the local (^-function and the Hadamard regularisation 
scheme in static spacetimes 

We intend to give a rigorous and well-known meaning to the formal definition (ITT]) in the 
Euclidean case relying upon an analytic continuation procedure, another procedure will 
be discussed later. We start by providing a definite meaning |Haw77] to the determinant 
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detP in the former identity ([T^ . If P is a tt, x n positive-definite Hermitian matrix witli 
eigenvalues < Ai < A2 < ■ ■ ■ < A„ the following elementary result is valid 



detP = e — ^1^=0 if C(s|P) := e C. (18) 

This trivial result can be generalized to non-negative self-adjoint operators whose spec- 
trum is discrete and each eigenspace has a finite dimension as it happens for the unique 
self-adjoint extension P of Laplace-Beltrami operators in compact Riemannian manifolds 
or perturbations of the form P = — A + in f llOp . Consider the series with s G C (the 
prime on the sum henceforth means that any possible null eigenvalues is omitted and, as 
is customary, we write P in place of P) 

as\P):=Y!x-\ (19) 
j 

As is known, the series converges absolutely to an analytic function if Re s is large enough. 
As in ( ITS]) , the idea |Haw77j is to define, once again, if the right-hand side exists: 

detP:=e-^l- and S,n := -'--^l., . (20) 

2 ds 

However, the function ( on the right-hand side is has to be understood as the analytic 
continuation of the function defined by the series (fT9|) in its convergence domain, since the 
series may diverge at s = - and this is the standard situation in the infinite-dimensional 
case! For P as in f IlOp the procedure works and (120]) makes sense (e.g. see [Mor99aj ). 

The idea can be further implemented giving rise to a (^-function procedure to compute 
the stress energy tensor, suitably interpreting f lTT]) . As a matter of fact |Mor98] . one takes 
the gfj,iy functional derivatives of the right hand side of f lT9|) : 

and, after an analytic continuation of Z^^{s,x\P) in a neighborhood of s = 0, define 

{■.T^u- {x)}n-= j-^ \s=o- (22) 

The arbitrary mass scale /z has to be introduced due to dimensional reasons. The outlined 
procedure works (the functional derivative can be computed giving a precise sense the the 
right-hand side of ( 1211) and the series converges to an analytic function for Re s sufficiently 
large) at least for compact Euclidean manifods with P as in (|TOl) as proved in |Mor98] . In 
particular it can be verified that (:T^jy: {x))q defined as in (1221) is conserved and produces 
the conformal anomaly. However the overall issue remains: how to relate (: T^^ : {x))q 
with the renormalized stress energy tensor in Lorentzian spacetime and for a state fl7 An 
answer was presented in |Mor98j and |Mor03j (see especially theorem 2.2 in the latter) 
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Theorem 6 Let (M, g) be a globally hyperbolic spacetime endowed with a global Killing 
time-like vector field normal to a compact spacelike Cauchy surface and a Klein-Gordon 
operator P in ^ where V does not depend on the Killing time. Consider a compact 
Euclidean section of the spacetime {ME,gE) obtained by: 

(a) a Wick analytic continuation with respect to the Killing time, and 

(b) an identification of the Euclidean time into Killing orbits of period /3 > 0. 

Let S be the unique Green function of the Euclidean Klein-Gordon operator associated 
to P by means of the Wick rotation when Pe is assumed to be strictly positive, and VL the 
Hadamard state associated to S . 

{'.Tfj^u- {^))n computed through the (-function procedure ^2^) coincides with the analytic 
continuation of Killing time in static coordinates (Wick rotation) of {-.T^^: {x))q computed 
through the point- splitting prescription as in ^ when the scales A and fi are chosen 
suitably. 



6 The well-posedness of the DeWitt-Schwinger regularisation prescription in 
Lorentzian spacetimes 

As described in subsection 13. 2[ one can formally obtain the expectation value of the reg- 
ularised stress-energy tensor by first regularising the effective action and then taking the 
functional derivative with respect to the metric. If such regularisation of the effective 
action leads to a diffeomorphism-invariant action, the associated stress-tensor expecta- 
tion value is automatically conserved |Wal84t app. E]. The original motivation behind 
the DeWitt-Schwinger point-splitting regularisation was indeed to proceed in this way, 
whereby the effective action is regularised by means of the formal asymptotic expansion 
of the heat kernel. However the actual computations in the seminal works |Chr76t IChr78j 
have been performed directly on the level of the stress-tensor, without the detour via the 
effective action. Nevertheless, we shall now repeat the formal steps of this detour before we 
proceed to give a reformulation of the regularisation prescription used in |Chr76i IChr78j 
which is both rigorous on general smooth Lorentzian spacetimes and takes into account 
the state dependence of (iT^;^:)^. 

To wit, inserting the heat kernal expansion (I14p in the definition of L^s fll3p . one finds 
that the contributions due to n = 0, n = 1, and n = 2 lead to divergent integrals with 
respect tosifx = y (or if x and y are lightlike related in the Lorentzian case). Particularly, 
we observe UV- divergences at the lower integration limit, while infrared divergences at the 

2 

upper limit are not present on account of ^. One therefore defines the renormalised 
effective Lagrangean in the Riemannian case as 

oo 

f 1 1 2 

Les, renix) := - ds -- — - c"™'" V a„(a;, a;)s" , (23) 
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and (:T^i/:)f7 by the functional derivative of the associated renormalised effective action 



S, 



eff,ren 



dqX Lgff 



M 



In the Lorentzian case, one would define these quantities by replacing s with is. We point 
out two important things. First, we know that the Hadamard coefficients f„ and, hence, 
a„ are covariant biscalars. Therefore, the renormalised effective Lagrangean is a scalar, 
and, consequently (:T^i/ :)q is automatically conserved. Secondly, it is clear that Lgg, ren(a;) 
has no state dependence whatsoever. This holds because we know that the coefficients a„ 
are completely specified by local curvature terms and m, ^. In fact, we have "lost" the 
state dependence by disregarding the smooth remainder term in the expansion of the 
heat kernel. Defining (: T^,^ :)n via Leg, ren(a;) therefore completely disregards the state 
dependence of (:T^y i)^. Apart from the appearance of the Hadamard coefficients, there 
does not seem to be a close relation to our definition of (: T^,y :)n in terms of applying a 
suitable bidifferential operator to the regularised two-point function and then taking 
the coinciding point limit. However, one can reformulate the above renormalisation of 
the effective action in the following way |Chr76irWa79j . One formally pulls the functional 
derivative with respect to the metric in the definition of {:T^^:)q under the integral with 
respect to s and then finds via additional formal steps 



6S, 



cff 



SS 



2 6^\detg\P{x,y) 



^/\detg\ 



6g 



/ii/ 



,/\d^\6g''''6(l)ix)5(l)iy) 



ds e 



-sP 



ST, 



flU 



5(l){x)6(l){y) 



ds e 



ds e 



-sP 



-sP 



= [P-'] {x,y) 

where the outer square brackets denote the coinciding point limit. In the above formal 
derivation, it has been used that the integral kernel P{x,y) = 6{x,y)Px of P is obtained 
as the second functional derivative of the classical action 5* with respect to the field and 
that the canonical differential operator D'^J^ we have considered in section 13.11 is nothing 
but the second functional derivative of the classical stress-energy tensor with respect to the 
field 0. In the context of renormalisation of the effective action in Lorentzian spacetimes, 
one usually considers P^^ to be the Feynman propagator A p (note that P~^ is not unique). 
Hence, the divergences of (: T^^ computed as above are interpreted to stem from the 
divergences of the Feynman propagator at coinciding points. To renormalise (: T^^, :)q 
in the Lorentzian case, one therefore inserts the DeWitt-Schwinger expansion (IMI) of the 
heat kernel in the integral expression for P~^ in terms of e~*^ to obtain 







^anix,y)s' 



n=0 
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n-1 



n=0 ^ / 

where the ^-prescription suitable for the Feyman propagator has been inserted and an 
integral identity for the modified Hankel function has been used. Expanding this in 
powers of cr, inserting = u, and removing the e-prescription from the regular terms, 
we find 

Ap(x,y) = — <^^— + log ^ ^ + + 



^■n'^ \a + ie ^\ 2 2 8 4 

rn?u bm^ua aia a2(T 
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and, inserting the relation between a„ and f„ as given in (fT7|) and ( [T6|) . we see explic- 
itly that, barring the different ^-prescription, lS.p{x^y) displays the Hadamard singularity 
structure. Again we point out that the 'correct' Feynman propagator is always state- 
dependent, while the above expression is manifestly state-independent, being essentially 
only a local curvature expression. Once more, this stems from the fact that one has 
disregarded the smooth (non-local) remainder in the expansion of the heat kernel. Note 
however, that, while this smooth remainder term is essentially well- understood in the 
Riemannian case, this does not hold in the Lorentzian there the whole DeWitt- 

Schwinger expansion is not rigorous. Hence, in contract to the Hadamard expansion of 
the two-point function of a Hadamard state, there does not seem to be a possibility to 
introduce the state-dependence in a rigorous way in the DeWitt-Schwinger renormalisa- 
tion as we have presented it up to now. Moreover, it is already visible from the few terms 
we have provided that, in term of the Hadamard series, the smooth term w of the above 
distribution contains inverse powers of the mass and, hence, diverges in the massless limit. 
Of course this is particularly also the case for the logarithmic terms, which displays the 
infrared singularity of the integrals with respect to s. 

By the well known distributional identities (see for instance |ReSi75t lFul89] ) 

1 1 

lim = V — h iirSix) , lim log(x + ie) = log |a;| -|- niQi—x) , 

where V denotes the principal value, one finds that [D'^^{x,y)[P~^]{x,y)] is a complex 
number, and one therefore has to consider its real-part as the 'correct' definition of (:T^i/ : 
)q. In terms of the Hadamard series we have discussed in subsection 13. H this corresponds 
to consider the symmetric part 

h'{x, y) := i {h{x, y) + h{y, x)) 

in the definition of the point-splitting prescription, where here and in the following we 
omit the emh^. Note that this encodes the same information as the full h in the coinciding 
point limit. 
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With the setup we have just described, the regularisation of {:T^u-)n goes as follows 
|Chr76] . Applying D'^^J^ to the real part of Ai;' given in terms of the DeWitt-Schwinger 
expansion, one finds that, as in the regularisation of the effective action, divergences come 
from the terms of order n = 0, n = 1, and n = 2. Hence, one identifies the divergent 
part of the stress-energy tensor as -D^™ applied to the first three terms in the expansion. 
However, these terms of course also contain smooth contributions, and in fact one has to 
take care to not subtract 'too much', otherwise one could spoil the covariant conservation 
of {:Tfj^i,:)n which was automatic in the renormalisation of the effective action. Namely, 
although in the latter renormalisation one has subtracted the first three terms of the 
DeWitt-Schwinger expansion as well, some of the subtractions have been zero on account 
of the vanishing of a in the coinciding point limit. As we now derive the series, we could 
accidentally introduce these vanishing terms since second derivatives of a do not vanish 
in the coinciding point limit present in Leg{x). However, as observed in |Chr76j . this can 
be avoided if one defines the divergent part of {:T^u-)n by applying to the real part 
of Af{x, y) and then discarding all terms which are proportional to inverse powers of the 
mass. Proceeding like this, Christensen has computed in |Chr76t IChr78] the divergent part 
of the quantum stress-energy tensor and has obtained the trace anomaly by computing 
the negative trace of the divergent part. This follows because the full expression of the 
stress-energy tensor must have vanishing trace in the conformally invariant case as, despite 
its divergence, it its completely given in terms of the real part of A^?, which in turn is a 
bisolution of the Klein-Gordon equation by its very construction in terms of the DeWitt- 
Schwinger series. Note that Christensen has introduced the massless limit by replacing 
the m? in the logarithmic divergence by an arbitrary scale A. This may seem rather ad- 
hoc, but from our point of view this is a reasonable procedure if we remember that we 
have defined the Hadamard series with an arbitrary length scale in the logarithm right 
from the start; moreover, as shown in |Wa79j . the scale can be consistently introduced as 
an "IR-regulator" . Finally, since the smooth terms proportional to inverse powers of the 
mass had been discarded to assure conservation, the massless limit could be performed in 
a meaningful way. 

Having reviewed the DeWitt-Schwinger point-splitting renormalisation of the stress- 
energy tensor, let us recapitulate its seeming disadvantages. 

a) It has been defined via an expansion of the heat kernel which is not well-defined in 
general curved Lorentzian spacetimes. 

b) It does not take into account the state dependence of {'■T^u:)q. 

c) It employs a Hadamard series whose smooth part w diverges in the massless limit. 

We shall now give a regularisation prescription which closely mimics the one of Chris- 
tensen, but disposes of the above three problems. 
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Theorem 7 Let VL2{x,y) := {(j){x)(j){y))n be the two point function of a Hadamard state 
Q, let 

^2i^, y) ■= 2 (^2(2;, y) + ^2(2/, x)) , 

and let us define for x and y in a geodesically convex neighbourhood 



-1- u 

hosix, y) ■= {'P- + loj 



6tt'^ I a 



\ \ 

2 2 8 4 4 



2 16 2 4 2m2 

:= ho{x,y) + h^{x,y) , 
1 0^2 

hm{x,y) := -^1^ , ho{x,y) := hDs{x,y) - hm{x,y) . 
Moreover, let us split the canonical bidifferential operator 

^KG,can ^ _ 3^)^.' V., - 2^ V. + ^G^. + {2^0, + (2^ - 1-) 4^'^,' " \m' 



as 

1 



Y\KG,can r-)0 1 y^m pim lim'^n piKG,can ptm 



The stress- energy tensor regularisation prescription defined as 

fulfils the Wald axioms of state-independence, local covariance and covariant conservation. 
Particularly, it displays the trace anomaly 

Proof. First of all, let us remark that the regularisation prescription is well-defined, as the 
relation between a„ and w„ given in (fTTl) and (fT6!) implies that — /ids is of class (the 
worst terms in ~ ^ds ^^re of the form cr^logo"). Additionally, the prescription fulfils 
the requirement of local covariance, since it only involves subtraction of objects given 
in terms of the DeWitt-Schwinger/Hadamard coefficients. Moreover, state- independence 
follows manifestly from the definition, as the modification of the canonical prescription is 
given in terms of D^j^hm, which is a state- independent term. 

To prove covariant conservation, we recall that, in the proof of theorem |H it has been 
implicitly computed that for any smooth biscalar B{x,y) the following relation holds 

[/}KG,can ^] ^ _ ^y^,p^ p] _ 
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Applying this to our current case, we find 

= - [v.. P. (fi^ - /ids) + V..P° hj\ = - [v.. {P.ni - P.h^s + P^hJ)] , 
where we have defined 

P° := -U + iR. 

A straightforward computation employing the Hadamard/DeWitt-Schwinger recursion 
relations yields 

Phus = P^hm ■ 

From this and the fact that Q2 naturally solves the Klein-Gordon equation, conservation 
follows. 

By Wald's results |Wal78a] . the above findings already imply that displays 
the "correct" trace anomaly. However, it is instructive to compute it explicitly. To 
this end, we obtain with steps similar to the ones already taken and using the implicit 
computational results obtained in the proof of theorem H] 

= - [(P, - m^) {n^^ - /ids) + Pl^h^] = [(]^ - /ids] = [fi^ - ho] - [a^] 

2vi - m vq-I — . 

□ 

A few comments on the result are in order. First, on practical grounds, the above result is 
really equivalent to the computation of Christensen in |Chr76l [Chr78j because the terms of 
the DeWitt-Schwinger expansion we have omitted are all proportional to o"^ and, hence, 
vanish upon application of the occurring differential operators in the coinciding point 
limit. In this sense, we have been able to put his results on firm grounds. Secondly, the 
modification term D^j^hm does not "simply cure the conservation anomaly". In this sense, 
the regularisation prescription just analysed differs from the one introduced in |Wal78a] 
and improved in |Mor03j in that it assures conservation in a different way. Namely, 
conservation does not follow by adding a term by hand or by modifying the classical 
stress energy tensor. In contrast, it follows by the explicit structure of the smooth term 
w in the DeWitt-Schwinger two-point function in combination with discarding specific 
terms proportional to inverse powers of the mass. Recall that the latter procedure has 
been motivated by analysing carefully the subtractions in the (non-rigorous) regularisation 
of the effective action. Finally, let us remark that the above prescription yields a smooth 
(: although Q2 ~ ^ds is only known to be twice-differentiable, because the non- 

smooth terms proportional to logo" all vanish in the coinciding point limit. 



[n'2 - ho] 
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7 Conclusions 



In this work we liave proven parts of Wald's conjecture in |Wa79j . namely, that the local 
C-function and the DeWitt-Schwinger regularisation prescriptions of the stress-energy 
tensor for quantum fields in curved spacetimes both satisfy the essential Wald axioms 
of state-independence, local covariance and covariant conservation, albeit only in static 
Lorentzian manifolds in the (^-function case. Moreover, we have shown how to reformulate 
the DeWitt-Schwinger prescription in such a way that it is able to take into account the 
full state dependence of the stress-energy tensor expectation value. Consequently, both 
prescriptions are physically meaningful for all we know and the results obtained by using 
them can be trusted although they do not seem to be rigorous at first glance. 

As a side note, it turns out that both the DeWitt-Schwinger and the Hadamard pre- 
scriptions require an equal amount of computation work in contrast to the seemingly 
widespread attitude that the latter is rigorous whereas the former is better suited for 
computations (see e.g. |DeFo06[ IDeFoOSj ) . Notwithstanding, the power of the Hadamard 
prescription is that can be directly generalised to obtain a regularisation of all field poly- 
nomials |BFK96l IBrFrOOl IHoWaOll IMorOSl IHoWaOSj ; this is done by means of the mathe- 
matic tools of "microfocal analysis" |Hor71t IDuHo72] IHor90] , which have been introduced 
to the framework of quantum field theory in curved spacetimes in |Rad96at IRad96b] . It 
is not clear how this can be done with the DeWitt-Schwinger prescription. 

Finally, we would like to mention a consistency problem of the semiclassical Einstein 
equation which is not solved in full generality by any regularisation prescription - this 
equation is only well-defined in general if one can associate to a each spacetime M a 
"unique" state Qm, because the spacetime is unknown prior to solving the semiclassical 
Einstein equation, but solving the equation is not possible before saying which Q is chosen 
for the evaluation of (:T^i/ :)n- However, it is known that such an assignment of a state to 
a spacetime is impossible |HoWa01] IBFV03t IFeVell] . Notwithstanding, this problem can 
be overcome if it is possible to restrict oneself to a limited class of spacetimes on which a 
unique state can be defined, this has been done in [PinlOj . 
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